Fast simulation of a quantum phase transition in an ion-trap realisable unitary map 
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We demonstrate a method of exploring the quantum critical point of the Ising universality class 
using unitary maps that have recently been demonstrated in ion trap quantum gates. We reverse the 
idea with which Feynman conceived quantum computing, and ask whether a realisable simulation 
corresponds to a physical system. We proceed to show that a specific simulation (a unitary map) is 
physically equivalent to a Hamiltonian that belongs to the same universality class as the transverse 
Ising Hamiltonian. We present experimental signatures, and numerical simulation for these in the 
six-qubit case. 



I. INTRODUCTION 

Feynman suggested that it is possible to simulate one 
quantum system with another^. However, we will turn 
this thesis around by posing the question of what sort 
of a system some unitary map on a quantum computer 
might correspond to. 

In particular, we examine the ion-trap model of quan- 
tum computing, and find that the unitary maps which 
have been realised on these correspond to the time evolu- 
tion of Hamiltonians which are linked closely to the Ising 
model. Finally, we consider the considerable theoretical 
body of work concerned with quantum phase transitions 
and renormalization group theory. This will later be the 
key to the problem of identifying a quantum phase tran- 
sition in a unitary map. 



A. Simulating Quantum Systems 

Feynman's first conception of quantum computingjpj] 
held the simulation of quantum systems as a key goal. 
Entanglement has been described as the quintessential 
feature of quantum mechanics 0. In general, the arbi- 
trary time evolution of a system is considered an NP- 
hardj^] problem, as memory and processing resources in- 
crease exponentially in the size of the problem, n, on a 
classical computer. It is only for extraordinarily simple 
systems, or ones for which there are strong symmetries, 
that such calculations are tractable. 

Feynman suggested that the problem could be reduced 
to one in polynomial time on a computer based on quan- 
tum principles. These include the ability of a quantum 
system to perform unitary operations on a set of quantum 
bits (qubits), and to exist in entangled states. Feynman 
showed that, in principle, it was possible to perform, in 
polynomial time, algorithms which were only possible in 
non-polynomial time on a classical computer. 

Since the original formulation of the problem, the ap- 
plication of quantum computing to classical problems 
has become more common. Several algorithms have 
been suggested, including the Deutch-Jozsa algorithm[4j, 
Shor's factorization algorithm^, and Grover's searching 



algorithm^. However, all of these systems are widely 
considered far removed from current experimental abili- 
ties. 

Recently, LloydQ revisited Feynman's original prob- 
lem, and showed that it was possible to implement the 
time evolution of an arbitrary spin Hamiltonian to a par- 
ticular precision, s, in polynomial time. The procedure 
essentially involves the decomposition of a Hamiltonian 
into realizable (local) unitary operations, and the time- 
wise stepping through a Hamiltonian to some arbitrary 
accuracy. 

It will be our desire to avoid such an abstracted simu- 
lation of a quantum system, and rather consider the pos- 
sibility of finding a quantum phase transition in a quan- 
tum algorithm naturally realizable with current quantum 
computing experimental hardware. In this way, we will 
essentially reverse the Feynman thesis, and conclude that 
quantum algorithms (or unitary maps) will correspond to 
the observables of some physical system. 



B. Ion Trap Quantum Computers and the Ising 
Model 



DiVincenzoy and Barenco et a/.^] has shown that 
single-site rotations and two-site controlled NOTs are 
universal for quantum computation. Further, the 
Sorensen-Molmer 10], phase gate ^lU]. and indeed almost 
any two-site entangling gate[l3 are universal. Hence, 
they will be able to affect any unitary transformation. 

Cirac and Zoller's paperjl^l on cold ion-trap quantum 
computers introduces the use of a spatially confined ion 
spin as a qubit, and the excitation of vibrational modes as 
a means of coupling qubits. Further, it has been shown 
that high fidelity state-preparation^3 and readout [iflj 
are feasible. 

Milburn has suggested a robust phase space scheme to 
use ion traps to simulate nonlinear interactions in spin 
systems^^l- A significant advantage of this scheme is 
that it does not require the cooling of vibrational states. 
The method involves the application of Raman pulses 
faster than the vibrational heating time, effectively de- 
coupling the effect of vibrational modes. In particular, 
Milburn shows that the evolution of a Hamiltonian of the 
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form 

H^nt = (1) 

may be achieved by a pulse sequence 

where X = and P = s-i^d expressions for 

and Kp given in Ref. [l7j . 

Further, Wineland's research group have recently 
demonstrated considerable success in achieving few- 
qubit interactions with this scheme. In particular, they 
present a two qubit phase gate, which has the form 

Ui>^Ui> , ITT) ^ ITT), 
UT)^e**UT) . ITi)^e'^ITi) 

which can be recast as 1^*) g^^x<^i^' <'i^'' |\i>^_ Apart 
from an uninteresting global additive phase, this may be 
considered to model the time evolution of a Hamiltonian 
of the form ai"'^a^'^^K 

Further, it is well known that single rotations in any 
basis, which correspond to the evolution of a spin opera- 
tor, are easily implementable on such an architecture [T^ . 
They result in unitary transformations of the form 

U single — 6 ^ (3) 

which can be implemented trivially through a single Ra- 
man pulse. 

Following Feynman's original intentions for quantum 
computing, one may consider the mapping of Hamilto- 
nian with such terms onto an ion-trap quantum com- 
puter. Turning this problem around, we will consider the 
properties of a unitary map composed of terms which can 
be experimentally implemented, and investigate their re- 
lationship with the transverse Ising spin chain. 

C. Quantum Phase Transitions and Universality 
Classes 

The quantum phase transition in the one dimensional 
transverse Ising model^^l is very well understood. The 
Hamiltonian is given by: 

N 

H,,»g-^Mi3(Ti") + Ja("V(") (4) 

n=l 

It is known that for an external field with interaction 
strength and local exchange interaction term with 



strength J, that a phase transition occurs for = ± J. 
One can intuitively consider the phase transition as a 
result of the incongruent symmetries between the two 
phases, which is reflected in the difference in behaviour 
of the two terms in the Hamiltonian under the transfor- 
mation cr„ — (T„. In the regime J > fiB, the system 

is in a ferromagnetic phase, with (^'^x^^'^ 0, and the 
system displays long range order. On the other hand, for 
J < ijlB, the system is paramagnetic, with^cri"^^ = 0, 

and there is no broken symmetry. 

Using arguments from renormalization group theory, 
we may place a great number of related problems into 
the same universality class and we may expect to see 
a similar phase transition occur in a number of related 
systems. 

II. THE MODEL 

In the following, we will put together the components 
introduced in Section ^ in a intuitive way. We consider 
the composition of the two unitary maps, similar to those 
demonstrated in Ref which corresponds to the com- 
position of the time evolution of two Hamiltonians. In 
form, it will look similar to the one-dimensional trans- 
verse Ising chain Hamiltonian. We will then apply the 
Jordan- Wigner transformation to this model to express 
the Hamiltonians in terms of non- interacting fermions. 
We are then able to perform a composition of operators 
in an SU (2) representation to yield a single Hamiltonian. 
We will flnd that this model is highly non-local. How- 
ever, using renormalization group theory concepts, it can 
be shown that the Hamiltonian belongs in the same uni- 
versality class as the transverse Ising chain. Hence, we 
conclude that our separated model has the same quan- 
tum phase transition as the transverse Ising chain, even 
though we have implemented the map in a much simpler 
way. 

A. Model Unitary Transformation and 
Experimental Realization 

It is natural to decompose the Ising Hamiltonian, 
Hi sing into two distinct parts: 

N 

ffx=xE'-i"^'-i"-''^ (5) 

N 

He^eY^ai-^ (6) 

n=l 

These parts are of even and odd symmetry under 
oVi — > —(^n, respectively. Unitary maps of the form 
|\I>) e'^^x.e have been realised experimentally. It is 
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impossible to perform them both at the same time with 
only single qubit rotations and two qubit gates, because 
they do not commute - the evolution of the combined 
Hamiltonian is not the composition of the evolutions of 



both Hamiltonians. Note however that terms a, 
do commute, and so 



(n](n+l) 



J (m) (m+1) 

and az cTz 



N 

n 



(")„(" + !) 



(7) 



— g ^j = l 3 ^ 



J - 



ale 



(13) 

they obey fermionic anti- 



It can be shown that 
commutation relations. 

We may understand these as an expression of domain 
wall creation and destruction. We can re-express U{xi^) 
with this new set of operators as 



is realisable in principle with current technology. 

The combined Hamiltonian may be approximated by 
Lloyd'sQ methods, which involves applying terms such 
as and ^Hq repeatedly, m times. However this 

requires a large overhead - instead we will consider the 
unitary map 



Finally, we will define the Fourier transformed versions 
of the fermion operators as 



Uix, 6) = e-'"^e-'"' = e''" ^ e-''^"^+"'^ (8) 

This map has been proposed by Milburn et al. 
as an easier unitary map to to simulate than the map 
which corresponds to the time evolution of transverse 
Ising chain Hamiltonian. We are interested in whether 
this mapping will have the same quantum phase transi- 
tion behaviour as the transverse Ising chain. 



= 7f E^'^e-'^ (16) 

- -^E^y"" (17) 

However, it is important to take note of the bound- 
ary terms. Strictly, in order to have Eqs. lEIl and Ijl2|l 
identical, we must make the identification|21j 



B. Jordan- Wigner Transformation 

We will follow Jordan and WignerpOj in using the fol- 
lowing definitions to introduce a new set of operators, a„, 
where 



An) 

X 

V 

Jn) 



1 - 2a„a]^ 



-i{a„ 



(9) 
(10) 

(11) 



where ui"'' , o-y"^ and di"'' take the form of the Pauli spin 
matrices in the |0),a|j |0) basis. From these definitions, 
the operators a„ and ojj can be shown to obey the fol- 
lowing relations: 



{ajj,a„}==l, al = 0, 



0, 



[al^, a„] = 0, [aj„, ajj = 0, [a„, a„] = 0,m^n 
With these definitions, our unitary map becomes 



tl2) 



We then introduce the following operators 



It may be argued that in the thermodynamic limit, this 
term will be irrelevant, and we may make the identifica- 
tion cjv+i — ci. 

Due to cyclic boundary conditions, we will require k 
to take the discrete values 



27rm 
k = — - — , m 



L L-2 



These operators satisfy fermion anti-commutation re- 
lations: 



{Ck,Cl} - Ski 

{Ck,Q} = {Clcjj^O 

Using the definitions of c„ and cjj , and the thermody- 
namic limit we can re-write U{x, d) as 



k-tsink(Clc\+CkC-k] 



18) 



where we require the thermodynamic limit so that the 
property 2^^ C^C^ = EJC^Cfe + Cl,C_fc) holds. 
To simplify matters, let us further define 
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Ak = x(2cos/cC|Cfc 
Bk = e{2ClCk-l) 
such that we may write 



isinA:(C|.Clfc 



CfcC_fc)Xl9) 
(20) 



Ug-^E.B. =nfcC/fc(x,0) (21) 



where Ukix,(^) 

We have now completely decoupled the problem, and 
may express the operators A/, and in the basis 
|0) , Cl |0) , Cl^ |0) , C^Clfe |0). It is be possible to find 
eigenstates of U (x, 0) in closed form in this basis. 



C. Combining 

However, it would be nice to be able to express U{x, 0) 
as a single exponential. While A]^ and Bk do not com- 
mute, it turns out that there is a faithful representation 
in SU{2), if we make the following definitions : 



(k) 



Hence, we can express Aj. = x(cosfc + a.^^) and B^ = 
(^Pk-i'k, where dk — x(sin fc, 0, cos fc) and Pk = ^(0,0,1). 



Sk where ei^r. 



(k) ,(k) ,,(fc) 



} have 



We have that [i^l^\i^ly:'^] ^ -2ie,, 
is the Levi-Civita symbol, so that {i^}' 
the same properties as the SU{2) matrices {cri, cr2, 0-3}. 
Relating the fermionic operators to SU{2) in this way 
was inspired by a similar approach in the theory of 
superconductors j^lj . 

SU (2) is closed under composition with a well under- 
stood composition relation, which we can now apply to 
our svstem[2^ 



be derived using quaternion composition[24- However, 
when using quaternions, special care has to be given to 
the double cover of 5*0(3) under SU{2). The second 
equality of Eq. Ij22|l defines an effective Hamiltonian, 
Hk, and we stress that Hk 7^ Ak + Bk because y4feand Bk 
do not commute. 

Hence, we have the final form of the decoupled, and 
combined transformation 



U{x,0)^nkUk{x,0)^e 



(26) 



Hence, we have found that the effective Hamiltonian , 
-ff, defined in Eq l|Sl is given by -ff = "^^^ Kkjkix^ 0)-^k 

We now check the limit x ~^ 0) which implies 
K ii^' {0, 0, cos fc sin 0}. Hence 



hand, 

Ik 



other 
and 



-^2ra^ay ^ ^e Y. k . Qu the 

in the limit ^ 0, k = jj^, 
= {sinx sin fc, 0, sinx cos fc}. Hence C/(x, 6*) — 
g-ixEfcSinfci.i-Hcosfci.3_ Thus, we retrieve the expected be- 
haviour in the limit as we turn off either the exchange or 
external field terms. 

Having expressed J7(x, d) in this form, it is now possi- 
ble to show that it directly corresponds to some physical 
Hamiltonian. We may perform a Bogoliubov transforma- 
tion by defining some fermion creation operator 



7fc7fc = lk{x,d)-^k 



(27) 



with associated energy, ek = Kk- Hence, we may consider 
our ground state as a vacuum state |0 >, and excitations 
as jI |0). It is important to note here that the excitations 
of lowest energy will occur at an extremum of Ck ■ We can 
show that this occurs at fc = 0, tt by noting that 



dck 
dk 



pKk _ duk drjk 
dk drjk dk 



(28) 



from which it follows that 



Uk{x.d) 
where 

Ikix.O) 

Kfc 



g-ixS-k -i^k ^-iSfik-t^k 



g-i cos fcg-iKfc7fc (x.S).^^ 



(22) 



(sin k cos 9 sin x, — sin k sin 6 sin x, (23) 
(sin 9 cos X + cos k cos 9 sin x) ) 



cos ^rjk 

cos 9 cos X ~ cos k sin 9 sin x 

k k 
cos^ — cob{9 + x) + sin^ — cos(6' — x) 



(24) 
(25) 



This composition has the simple physical interpreta- 
tion of two rotations being composed, and the result can 



drjk 
dk 



sinfcsin6'sinx|j,^o TT ~ ^ (29) 



fe=0,7r 



Hence, the elementary excitations will be for fc = or 
k = TT, whichever corresponds to a lower energy. 



D. Closed-form Hamiltonian 

We can now work backwards from our expression for 
U{x,9) to a single combined Hamiltonian. The results 
here will only be valid in the thermodynamic limit, which 
we have assumed in the previous section. Before doing 
so, we should present a list of identities which will prove 
to be useful. 
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mfc (fe) ^ jS^^JJ r ^ 

k n 
k n 



^ sin(afc)z^f ^ = ^ 4^+^ - c„c„+a 
^sin(afc)i/^'^^ = i^(4cj,+^ +c„c„+a) 



Now, recall that 



(30) 



However, Kk is an even function of k, and so can be 
expanded in terms of a Fourier series involving cos k. Let 
us write 



Kk = ai cos(Zfc) 
1=0 



(31) 



^^^^ _ ^-iJ2k,l °-l cos(ifc)7fc(x,e).i7fc ^g2) 

We will now substitute our expression for 7^ and ex- 
pand. 



^ g-i(Ai+A2+A3) 



where 



Ai — cos ^ sin X a; sin k cos{lk)v'^^ 

k,l 

A2 = — sin sin x a; sin fc cos{lk)v^^ 

kd 

A3 = sin6'cosx^ a; cos(Zfc)i^3'^^ 

+ cos 9 sin X a; cos fc cos{lk)v^^ 



(33) 
(34) 
(35) 



(36) 
(37) 
(38) 



k,l 



and the sum over / ranges 1,2,3. 
We may rewrite this as 



Ai = coses\nx^ai[s\i\{l + l)k-sin{l-l)k]v''^\m) 

k,l 

= cos6'sinx[ao(cjj4i+i - c„c„+i) (40) 
A2 = — sin 9 sin x ^ ai[sin(? + l)fc - sin(/ - l)fc];/Pl) 

k,l 

= -«sin6'sinx^[ao(cJjcJj^;^ + c„c„+i) (42) 
sr^ (a/+i - a;_i) , + t , m 

+ 2^ 2 ^^nCn+l + C„C„+i)J 

A3 — sin 6> cos X a; cos{lk)iy^^'''' 



k,l 



- COS 9 sinx ai [cos(Z + 1)A: + cos(Z — l)fc] 



k.l 



(43) 



,(fc) 



(44) 



sin6'cosx^[a/(c|'jC„+; - c„4_^;)] + 
I 

+ cos6'sinx^[ao(cj,c„+i - Cnc\^-^) 

n,l 

(a;+i-a;_i)) ^ t M 

I 2 ISi'-n+i '-nSi+iJJ 



Therefore, the quantum spin chain Hamiltonian H 
which represents a physical system corresponding to the 
separated unitary map l|Hl is highly non-local. Terms 

such as c„c]j^Q — a„e^*'^^J=" °'^'^^o,n+a for a > 1 will 
not only involve a„ , a„+a , but also Cm and c\^yn < m < 
n + a. 

If we define 



An) 



(") , ■ (») 



(n) . (n) 
(„) _ nl ~i<^y 



we can write c„c„_,_^ as cr^ e ^3=^ ^ (t_ , ex- 
plicitly showing the dependence on non-neighbouring 
spins. 



E. Range of the Interactions 

To be in the universality class of the Ising model, we 
would expect the non-local terms, ai to decrease expo- 
nentially with separation I. Thus, when viewed at larger 
length scales, the non-local terms would become irrel- 
evant. The behaviour of a„ for a variety of = x is 
calculated numerically and presented in Fig. ^ The case 
for 6* ^ X is similar, and displays the same exponential 
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Figure 1: The behaviour of Fourier coeflEicients defined 
in Eq. E5| for n = 1, 2. ..10 for a variety of 6 and X- Note that 



the larger X ~ f and \9 - 



the larger the decay. The ex- 



ponential decay in these coefficients implies that the interac- 
tions in the Hamiltonian are short-ranged, and suggests that 
renormalization techniques should be highly effective in this 
model. 



decrease in a; with I. Hence, we may naturally expect 
nearest-neighbour interactions to be the most important 
interactions in this model - an idea which we will make 
concrete in the following section. 

Deriving an analytic expression for a„ seems to be very 
difficult. However it is possible to show a general depen- 
dence on 9 and x of the form sin" 9 sin" x for Par- 
ticular n. We can expand Kk in terms of Jjk, where ijk is 
defined in Eq. as 



p=0 ' p=0 



(45) 



where F is the Euler gamma function. In turn, if can be 
expressed as a series in terms of cos'' k as 



r/P = ^ ("^"j (cos0cosx)''"''(-cosfcsin6lsinx)(^6) 

P 

= dp, 5 cos* k 



q=Q 



Finally, we can express cos' k in terms of cosrfc. For 
the case of q even, this is 



cos*.= ± + c (47) 

9 



Cq^r COSrfc 



r=0 



where c is an unenlightening constant, and a similar ex- 
pression holds for p — q odd. 



We can combine Equations . and ijTzjl to yield 
an expression for Kk 



oo oo p g 

Kfe = X! ^P^^ = X! X! X! ^p'^p-.i^q.-r <^os rfc (48) 

p=0 p=0 g=0 r=0 

from which we can read the coefficient a/ of cos Ik as 



oc oo p 

p=0 p=0 g=0 



(49) 



However, Cg,/ is only non-zero for q > I, and so we can 
replace Y^^^q with J2q=i to yield 



oo p 



q,l 



(50) 



p=0 



p=0 q=i 



which involves terms in dp^q oc (sin0sinx)*(cos6'cosx)* 
for q > I. Hence, for any given ai, there is a behaviour 
proportional to (sin 6* sin x)'. 

Of physical interest is the behaviour of ai with respect 
to I for a given 9 and x- We will attempt to factor out 
any behaviour in / by noting that since dpq,eq,i < 1, we 
can write 



(51) 



ai = Cp dp^qCqJ 

p=0 q=l 
oo p p 

- X! XI ^P'? X! ^9^' 
p=0 q=l q=l 

However, if we turn our attention to the sum over dp 
we can construct a further limit on ai 



J2dp,q = ) (cos6lcosx)^"'^(-cosfcsin6'sinx)(^2) 

q=l q=l \ / 

^ E ( ,^;) (cos0cosx)''-''-'(sin0sinx)(^*j 

P~l / _ I \ 

- ' j (cos6lcosx)P"'"'^'(sin6isinx)'('5^4) 

< (sin 61 sin x)' (cos 6* cos X + sin sin x)''"' (55) 

< (sin 61 sin x)' (56) 

Leading to our strictest inequality for a/, showing an 
exponential decay in I 



oo p 

ai < (sin 9 sin x)' Cp Sg.i 

p=0 q=l 



(57) 



< (sin6'sinx)'(/^ Cppmaxcg./) (58) 



p=0 
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The case 9,x ^ j has asymptotically constant a; . 
However, for 0, x 7^ f j we can say that 



' sm{0 + X) ox ox 



ai<ke ^ ^0 as ^ oo (59) -2(6* + 

where /i — ln(sin^? sinx)- Thus the terms in our model jg some constant 

H has only short range interactions. Applying the transformation ^ e^i*, we have that 



F. Renormalizing the Hamiltonian 

We have seen that the Hamiltonian, H, associated 
with the complete unitary U is very complicated and in- 
volves non-local interactions. Hence, we are forced to use 
renormalization group methods to extract the interesting 
physics from this case. 

Consider the continuum limit of the Hamiltonian, H , 
which is applicable in the thermodynamic limit. Near 
criticality, the physics will be driven by long-wavelength 
effects, which suggest that the wavevector k will be small. 
The relevant excitations at low temperature happen at 
the extremum of e^, which we have shown occurs at 
k — Q. Hence, in the continuum limit, we can consider 
only low lying states, near k — Q. Under these approxi- 
mations, our Hamiltonian becomes 



(60) 



7fe = (fc cos6' sinx, — fc sin6' sinx, sin(0 -I- x)^61) 
+ X 



sin(6l + x) 

which yields in terms if the fermion operators 



(62) 



H = En 



dx[ 



(6* -Fx) sin X. id'i'^ 9*,, , 
sm{0 + x) ^ dx dx'^ ' 



-2(0 + X)*^^] 

If we do not perform this transformation, we will have 
terms of the form / dx^^ + "if^ in the Hamiltonian. 
One can show that these terms correspond to interac- 
tions of the form cri"''^^"^^'' -I- ay"''' ai"^^\ Chapter 4 of 
Ref. [2^ has further details regarding these sorts of chiral 
symmetries in systems. 

One can show that the Lagrangian corresponding to 
this Hamiltonian will then be 



dr 
-2(0 4 



+ X) sinXj^^t.^*^ 



sin 

t 



id + x) 



dx 



where r is imaginary time. 

Now we introduce the crucial step where we consider- 
ing the effect of scaling the problem. If we consider the 
effect of viewing the problem at a scale more coarse in 
space, and 5^' more coarse in time, we can introduce the 
new variables 



+2(9 + x)ClCk 
Defining the continuum Fermi field[2^ as 



^'(cci 



(64) 



where a is the lattice spacing, ^'(a;) satisfies the usual 
anti-commutation relation {'^'{x),'^^ {x')} — 6{x — x'). 
Note that we can replace the sum over k with an integral, 
by making the substitution 



dx 



Further, we expand the terms C|clj. and CkC-k into 

terms of first order gradients and ^ through the 
use of identities (Chapter 4 of Ref. ^3|)- This yields 



x' 


= x5-^ 


1 

T 









(68) 
(69) 
(70) 



We choose the value of the dynamic critical expo- 
nent, z, to be identically equal to 1, corresponding to 
an isotropy between space and time, in order to leave the 
velocity-like coefficients of ^'^^r and unchanged. 

For criticality to hold, these scaling conditions must 
leave the Lagrangian unchanged. This occurs only when 
the quantity -I- x is identically zero. This happens for 
9 = —X, exactly as in the transverse Ising model. 

Formally, if we write 



97 



u(\E' 



dx 



) + A*t^ (71) 

ox 



we will require that 



A' 

u' 



A6' 
u 



(72) 
(73) 
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implying that the scahng dimension of the term 
A\l/'f\l>,dim(A) = 1. Hence, u, and A are relevant pa- 
rameters, having non-negative scahng factors. 

If we include second (or higher) order effects in k, we 
will include terms of the form A'^-t^ or A"*t^||^ 
(or higher derivatives). From a simple analysis, one can 
show that the parameters A' and A" are irrelevant, as 
the scahng dimensions a.rejl^ 

dim(A') = -1 , dini(A") = -2 

Recall the Ising chain in a transverse field, Eq. 1^, 
given by 

JV 

HIs^ng = A^B^") + Ja^") ^^"+1) (74) 

n=l 

= J2'2BClCk + 2J COS kClCk (75) 

k 

+iJsmk{Clcl^ + CkC-k) 
The Lagrangian for this model has the form 



A. Ground State Energy 

Recall that we have the following unitary map which 
describes the system in terms of non-interacting fermions: 

By inspection, the corresponding Hamiltonian is given 
by: 

= ^Kfe7fe(x,0).j/fc (80) 
fe 

We note that eigenstates of H are simply products of 
the eigenstates of H^, where Hk = '«fe7fe(Xj ^)-^- In the 
basis |0) , Cl |0) , Clfe |0) , Cld^ |0), there will be four 

complex eigenvalues — 1,2,3,4}, with arguments 

i = 1, 2, 3, 4}. Let us define the argument of a total 

system state, fi*^*' = J2k^k ^^ where we form an eigen- 
state of H from matching eigenstates of Hk ■ Physically, 
we associate the argument of this state with energy. 
Now, we will consider the mapping 



s^na = ^'t|^+2(B-|- J)*t^- J(*t 



dx Ox 



Hence, we may make an association between the two 
models through the mapping 



J 



[e + x) sinx 



B 



s\n{9 + x) 

J = e + x 



(77) 
(78) 



Thus we conclude that our continuum Hamiltonian H 
belongs in the same universality class as the transverse 
Ising Hamiltonian, Hi sing- Hence it is possible to access 
the physical properties at criticality of this well known 
model in a very straight forward manner. The crucial as- 
sumptions which have been made are the assumption of 
operation in the thermodynamic limit, and the low tem- 
perature (and hence small k excitation) regime, both of 
which are necessary for renormalization to work. We will 
show in the next section that for moderate N ^ the signa- 
tures of quantum phase transitions are still observable. 



X ^ f cos (p , U r sm 

where (j) can be considered as a relative strength between 
exchange and field coupling terms, and r is an overall 
strength. The Ising criticality condition 9 — ±x is now 

In Fig. El we observe an sharp peak in the second 
derivative of the ground state energy, fj'^', with respect 
to (j). In the thermodynamic limit, this would become 
a singularity, indicating a second order phase transition. 
Further, we observe that this condition occurs for = 
±j,±^, which corresponds to the Ising transition. 

We will demonstrate the nature of this singularity ex- 
plicitly, by noting that the argument of the ground state 
energy is 



j^^-* = — (cos k cos sin x + sin cos x)^ + (sin k sin(^J.3 



cos rjk = E, 



k 



(82) 



The next two eigenvalues are equal to 1, and hence 
^k ~ ^k ~ ^- highest excited eigenstate has 



,(4) _ 



,(1) 



by symmetry. 



Substituting ^ r cos (f) and x ^ r sin </>, we evaluate 



III. SIGNATURES OF A QUANTUM PHASE 
TRANSITION 



The first experimental realisations of a quantum sim- 
ulations will perhaps be seen on ion trap quantum com- 
puters. In this section we review a few basic experimental 
signatures which may be seen in an ion trap laboratory. 



2N 



(cos k ~ l)dk 



(cos^ — 



72 



cos k siv? 



We find that the residue of the integrand is 



ir/V2' 



and hence has a -r singularity. We can now conclude that 



9 



— 1^ — 



in the limit as A'' — > oo, the value of 

will be infinite, with a logarithmic singularity. 

Alternatively, following Ref. and expressing Ek as 



Ek = cos ^ rjk 



(84) 



k k 
= cos^^(cos^ — cos{6 + x) + sin^ — cos(6' — x)) 

Bunder and McKenziefjll note that for some wave vec- 
tor fc, Ek — 0, corresponding to a vanishing energy gap in 
the system. When Ek is expressed as Eq Ij84|l . it is clear 
that there will be no energy gap for fc = if 6* = — x and 
for fc = TT if 6* = X- Without loss of generality, we can 
consider the k — the other is symmetric. Since 

the relevant excitations are at fc ~ 0, we may use Eq Ij84|l 
to expand El as a series in k 



Et (O + xr + k 

= e^ + fc^C^ 



2(^^ + x)(cos( 



X - COS ( 



2 sin(6' + x) 
and the ground state will have energy 

Eo ^:^j\k^. 



■85) 



(86) 



where kc is a cutoff wavevector. While analytical solu- 
tions are possible, it is unenHghtening to solve this prob- 
lem. Instead, we can set out to determine the behaviour 
of the energy with respect to a variable by using the 
indefinite integral 



E„ = - I ae^Eo 



(87) 



Carrying this out we obtain 



Eo = J 

-1 
47r 
-1 



,_d_ 

de 

dC 



:Eq 



Kc 1 



Thus we confirm the logarithmic nature of this sin- 
gularity, and find that Eg _^2-a ^j^gj-g .^^g have the 
value of the critical exponent = 0+. This is the same 
behaviour as that found in the transverse Ising modelfj^]. 
which we expect by their inclusion in the same universal- 
ity class. 

One can see the behaviour of ' with respect to (f) 
in Fig. |2lfor TV = 200 and r = 1.9. There exists a quan 
tum phase transition at 4 
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Figure 2: The second derivative of the phase of the eigenvalues 
for the model Hamiltonian H (solid), based on Eq. 12211 . and 
for the transverse Ising Hamiltonian, Hiaing (dashed), as a 
function oi <j) — tan~^ ^. Note that both show singularities 
&t 6 = ±x- 1^^ ~ '^^^ chosen to be 1.9 so as to 

highlight the differences between the plots. For smaller 6 and 
X, the commutator between Hg and becomes small, and 
the model becomes asymptotically closer to the Ising model. 
Fig. Ushows that the smaller 9 and x, the faster the model's 
non-neighbour parameters decay. 
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Figure 3: The second derivative of the ground state energy of 
a 6 qubit model as a function of = tan~^ |-, for r = 1.9. The 
maximum value is attained at = ^, the value at which a 
quantum phase transition occurs in the thermodynamic limit. 
While we see a strong maximum, in the thermodynamic limit, 
we expect to see a singularity. 




the singularity in 



This numerical modeling cor- 



±|-,±^ as evidenced by 



responds to our theoretical expectation for the positions 
of the phase transitions. For a finite set of qubits, one 
can clearly see the peak in the second derivative of the 
energy with respect to 4> in Fig. 



B. Entanglement 

It has recently been shown that entanglement scales 
near a quantum critical point (2^ l2^. Quantum phase 
transitions are driven by quantum fiuctuations.Igj, and 
entanglement is a natural manner for non-local effects 
to manifest themselves. As entanglement is a physical 
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Figure 4: The derivative of the nearest neighbour entangle- 
ment, with respect to ^ for a 6 qubit model. Note that this has 
a maximum very close to the critical point (p = j. We expect 
to see the nearest neighbour concurrence vary as log \ (t> — 
in the thermodynamic limit. 



resource, it may be directly measured, by a number of 
schemes |29ll3nl |. 

We denote the nearest neighbour entanglement in the 
ground state hy £. In the transverse Ising model, we see 
that the derivative of entanglement with respect to 0, 
near criticality to scale as a function of {(f) — ipd- Since 
we are in the same universality class, we expect to see 
identical behaviour in this model. However, experimen- 
tally, isolating the ground state of the system to observe 
this may be very difficult. 



C. Spectroscopic Measurement of Eigenvalues 

Experimentally, it is possible that an ion trap may be 
used to implement the unitary map of the form 

followed by another map of the form 

Let us introduce the notation 



where [vj/f™)^ is the state after we repeat this unitary 
map, U, m times. 

The state l^*) can be decomposed as 



(89) 
(90) 



where |0„) are the eigenstates of U , with "energy" En. 
Without loss of generality, we assume these are ordered 
with 



< E, 



where M = 2^. 

We can proceed to measure U"^ |^) in some set of basis 
states, \i), which will typically be binary computational 
basis states. Hence, we can measure 



imErj 



(91) 



(92) 



If we perform a Fourier transform of |(?|C/™|4')| over 
TO, we expect to see peaks around the allowable transi- 
tion energies £"„ — En' . A numerical simulation of this is 
shown in Fig. El for 4 qubits. Let us define 



En = 



m=0 



(93) 



to be these Fourier components. 

Since we are considering Unitary maps, and not Hamil- 
tonians, we can only determine the eigenvalues of U to 
within an additive constant of 2tt. Hence, in order that 
the Fourier components are not aliased (that the energy 
levels do not "wrap around" on themselves), we require 
the ground state to have an energy Eq > — tt, and the 
highest excited state to have an energy E2N < n. Since 
the energy is a function which scales with 0{N, 0, x), we 
require the condition max(|f?| , |x|) < -^j^ where kmt is 
0(1). Keeping 9 and x small in this manner will en- 
sure that the energies will be resolvable uniquely by the 
Fourier transform. 

If we have a given set of energy eigenvalues, 
{£^0)^1 • ■ -Em-i}, we can form the set of energy dif- 
ferences, {Eij = Ei — Ej^ . We can then calculate the 
Fourier transform of these differences, and compare our 
measured spectrum with the calculated spectrum. If we 
have n 3> M, then the problem is over determined, and 
we can apply a least-squares method to reconstruct the 
original energy spectrum (to within an additive constant, 
and global sign chan ge). We may apply the Levenberg 
Marquardt algorithm|3lL l^^ . to perform this recon- 
struction in polynomial time with an initial guess at the 
set of energy eigenvalues. 

One could change the value of |- over many experi- 
ments to tune the system through the critical coupling. 
In the thermodynamic limit, the energy gap to the first 
excited state would vanish at criticality, but we see in 
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Figure 5: Several of the lowest energy eigenvalues for the 6 
qubit model. One can clearly see the excitation gap closing 
as <f> approaches ^. In the thermodynamic limit, we expect 
the gap to be identically zero at = . However, we can still 
observe the energy gap, A, behaving as A ~ |(^ — with 
some higher order corrections. 



the relation 

A~|0-(/.,r (94) 

with 7 = 1, as we expect from a standard treatment of 
the transverse Ising problem jT3|. 

1. Controlled- U Spectroscopy 

The above method requires knowledge of the approxi- 
mate values of | {i\U"^\'^) \ , which means that a measure- 
ment with result \i) must be achieved multiple times. 
In work by Miquel et al^^, it has been shown that 
spectroscopy can be achieved much more easily by im- 
plementing a controlled-U gate, and measuring only a 
single qubit We can achieve this by using an ancil- 
lary qubit, and express the controlled-U operation as 



Fourier 
Component , m 




Figure 6: The Fourier transform of |(^|^7'"|*>|^Kn,2048, is 
shown as a vertical density, as function of the horizontal co- 
ordinate 9, for a fixed x ~ 0-2, and 4 qubits. 4 qubits are 
chosen so as to provide a complex, but not confusing diagram. 
The white bands indicate a large Fourier component. The 
superimposed grey lines show all energy difi^erences - note that 
some of these are disallowed. In this case, 2048 samples are 
used in the Fourier series, and simulations are taken in steps 
of 0.01 in 6. From this diagram, we can see the energy gap 
between the ground state and first excited state approaching 
zero. We can also see a level crossing, where one of the grey 
lines is reflected through the origin at S ~ 0.25. 



Fig. I^that the condition is not strictly met for a finite 
number of qubits. 

One can use this to show that the energy gap. A, for 
the excitation from the ground to first excited state obeys 



CU : {i\ ® (*| ^ {i\^U'{<f\ 

where {i\ can be either (0|, which takes to (^'|, or (1|, 
which takes {^\ to 

If we do a weak measurement on the control bit, we 
yield the result 

{<j,)^n[Tr{Up)] , K) = 3 [rr(C/p)] 

where p is the density matrix corresponding to the state 
has been prepared in. If we prepare it in the mixed state 
given by p = 1/2" where / is the identity operator, we 
yield (ctz) — [Tr{U)] /N, which is proportional to the 
sum of the eigenvalues oiU. If we repeat this for U"^ for a 
variety of m, we can use the method above to reconstruct 
the energy level diagram. 

Further, Miquel et a/. [34(1 propose a scheme using the 
quantum Fourier transform to probe specific regions of 
the spectrum of the eigenvalues of U. This is achieved by 
introducing an effective time scale into U, and exploiting 
the conjugacy of energy and time. 

D. Phase Estimation Algorithm 

In work done by Abrams and Lloyd|36i, ^^'^ f^^" 
ther explored in an ion trap context by Travaglione and 
Milburn^3|) it has been shown that it is possible to esti- 
mate the eigenvalues associated with any unitary trans- 
formation, U. These correspond directly to the energy 
eigenvalues of the equivalent Hamiltonian, which we are 
interested in. The scheme also yields an approximate 
eigenvector with high probability. 

Starting with a mixed index state \j)i, and the state 
of the target system j^*), we perform the transformation 

A{U) : |j)/|*>T^ |j)/®C/^'|*)t 
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followed by a Fourier transformation on the index reg- 
ister. Measuring the index register will then yield, with 
high probability, an approximate eigenvector of U in the 
target state, and information about the phase of the 
eigenvalue of U in the index register. 

Note, however, that use is made of an index register, 
which is at least the same size as the system of interest. 
This makes it a much more difficult problem to conquer 
experimentally, as a system twice as big will be much 
more prone to decoherence. In ion trap implementations, 
trapping twice as many ions will also be more difficult. 
While this technique is superior to the spectroscopic mea- 
surements suggested in the previous section, scalability 
issues may keep it from being experimentally feasible for 
some time. 



IV. CONCLUSION 

We have presented a number of key ideas which will 
drive our search for a quantum phase transition in a 
system which is implementable on an ion-trap quantum 
computer in a natural way. 

We have taken the Feynman thesis and turned it 
around, to ask what might happen if we have some 



implementable unitary transformation. The two fields 
which will help to answer this question have been intro- 
duced - namely, the ion-trap quantum architecture which 
will provide our unitary transformation, and the tools 
of renormalization group theory. In this paper, we have 
shown that the Hamiltonian corresponding to a separated 
Ising map belongs to the same universality class as the 
transverse Ising model. Further, the map presented here 
is realisable in a very natural way on an ion-trap quan- 
tum computing architecture. 

We have also suggested some experimental signatures 
, including ground state energy and entanglement, and 
spectroscopic information which may lead to the recon- 
struction of the energy spectrum. 

After completion of this work, we became aware of 
some other work on simulating quantum phase transi- 
tions in ion traps ^38j 
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